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Abstract 

A rapidly rotating Bose-Einstein condensate in a symmetric two-dimensional harmonic trap 
can be described with the lowest Landau-level set of single-particle states. The condensate wave 
function i/j{x, y) is a Gaussian oc exp(— r^/2), multiplied by an analytic function f{z) of the complex 
variable z = x + iy. The criterion for a quantum phase transition to a non-superfluid correlated 
many-body state is usually expressed in terms of the ratio of the number of particles to the 
number of vortices. Here, a similar description applies to a rapidly rotating non-symmetric two- 
dimensional trap with arbitrary quadratic anisotropy (a;^ < uj"^). The corresponding condensate 
wave function ip{x, y) is a complex anisotropic Gaussian with a phase proportional to xy, multiplied 
by an analytic function f{z), where z = x+i[3-y is a stretched complex variable and < /3_ < 1 is a 
real parameter that depends on the trap anisotropy and the rotation frequency. Both in the mean- 
field Thomas-Fermi approximation and in the mean-field lowest Landau level approximation with 
many visible vortices, an anisotropic parabolic density profile minimizes the energy. An elongated 
condensate grows along the soft trap direction yet ultimately shrinks along the tight trap direction. 
The criterion for the quantum phase transition to a correlated state is generalized (1) in terms of 
N/ Lz , which suggests that a non-symmetric trap should make it easier to observe this transition or 
(2) in terms of a "fragmented" correlated state, which suggests that a non-symmetric trap should 
make it harder to observe this transition. An alternative scenario involves a crossover to a quasi 
one-dimensional condensate without visible vortices, as suggested by Aftalion et al, Phys. Rev. A 
79, 011603(R) (2009). 

PACS numbers: 03.75.Hh, 05.30.Jp, 67.85.De 
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I. INTRODUCTION 



Most theoretical studies of rapidly rotating Bose-Einstein condensates (BECs) assume a 
symmetric harmonic trap with = ujy = >jj]_ In this C3iSG, clS the rotation 

frequency f2 approaches the trap frequency tu^, the iV-particle condensate expands radially 
to a large radius i?o, leading to a reduced central density. Quantized vortices provide the 
angular momentum with an essentially uniform vortex density ^ T{h/{MVL), where M 
is the particle mass. 

The situation can be quite different for a non-symmetric harmonic trap with uJl. < u^, 

nnnn ^ ^ , 

13, la, la, llOj. In this case, the condensate expands along the weak trap direction (here x), 



which enhances the moment of inertia and hence induces extra angular momentum. The 
energy functional in the rotating frame involves E — flL^, so that the optimal ground state 
for fixed Q tends to maximize L^. This description emphasizes the role of the deformation 
of the condensate. In some situations, the sequence of states with increasing Q can even 
reduce the number of vortices in favor of increased elongation (sj. 

In Sec. II, I briefly summarize the situation for a rapidly rotating symmetric trap, when 
the lowest Landau level (LLL) forms a nearly degenerate set of noninteracting single-particle 
states. These states serve as a variational basis for a condensate wave function i/jlll js,!^!. 
The criterion for a quantum phase transition (QPT) to a correlated non-superfluid many- 
body state is usually expressed as the "fllling fraction" N/N„, where N is the total number 

n n nri n n 

of particles and Ny is the number of vortices [HI, 112|, 113|, IM, 115|, 116[. Section III then 



discusses the behavior in a non-symmetric rotating trap and relies on the corresponding but 
more complicated LLL single-particle states in a rotating non-symmetric trap j^, |lo|. The 
rapid rotation induces an expansion along the weak trap direction and a contraction along 
the strong trap direction, increasing the moment of inertia and the angular momentum L^. 
Section IV uses these results to estimate the critical angular velocity for the quantum phase 
transition (QPT) in a non-symmetric trap. One proposed approach relies on the ratio of the 
total number to the angular momentum per particle Lz and suggests that it may be easier 
to achieve this QPT in such a non-symmetric system. A second proposed approach relies on a 
comparison with a "fragmented" correlated many-particle state 17| and suggests that it r nay 



be harder to achieve this QPT in such a non-symmetric system. A still different scenario 



id 



suggests an effectively one-dimensional condensate with a fixed minimum dimension in the 
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strong trap direction and few or no visible vortices (see also 18|, ll9|). Section V considers 
this possible crossover in some detail. 



II. REVIEW OF QUANTUM PHASE TRANSITION FOR SYMMETRIC TRAP 

Consider a dilute Bose-Einstein condensate with atoms, each with mass M, in a sym- 
metric two-dimensional harmonic trap with V{r) = ^Muj^r"^ = |Mu;^(x^ + ?/^). Apply a ro- 
tation in the positive sense with angular velocity Q. In the limit of rapid rotation {Q —>■ uj±), 
the radial expansion renders the condensate effectively two dimensional, and it is convenient 
to assume that the condensate is uniform in the z direction with thickness Z. The original 
three-dimensional condensate wave function can be rescaled as "^{r, z) = \J N/ Z iljir), where 
ip obeys the normalization condition j (Pr |?/'(r)p = 1. I use dimensionless variables with 
uj± and >yti/{Muj±) as units of frequency and length. In this way, the one-body Hamiltonian 
in the rotating frame is 

H', = ^ + --QL^, (1) 

where = z ■ r X p = xpy — ypx is the z component of angular momentum and p = 
— iV. With this choice of normalization, the Gross-Pitaevskii energy functional per particle 
becomes 



where g2d = 4:7ras/Z is the two-dimensional coupling constant with the s-wave scattering 
length (assumed positive) [2^. 



A. Mean-field Thomas-Fermi regime 

The mean-field Thomas-Fermi regime occurs when the rotating condensate has a dense 
vortex array with well-separated vortex cores. In this case, the density variations can be 
neglected, and the dimensionless momentum operator becomes the superfluid velocity v, 
which closely approximates solid-body rotation = fl X r. Minimization of the resulting 
approximate Eq. ([2]) with respect to the density yields the Thomas-Fermi (TF) profile 

mr)\' = n{0)(l-^), (3) 



3 



where n(0) = fi/{Ng2d) = 2/(7ri?y^) is the central density, /i is the chemical potential that 
enforces the normalization condition, and the condensate radius is given by 



As expected physically, the condensate expands with increasing rotation, and the central 
density correspondingly decreases. The criterion for the validity of the mean-field TF ap- 
proximation is fi = n{0)Ng2d ~ 1, which here implies 



^{1-^')>1. (5) 



B. Mean-field lowest Landau level regime 

In the opposite limit n{0)Ng2d ~ 1, the system enters the mean-field lowest Landau level 
(LLL) regime, when the vortex cores start to overlap. The density variation now becomes 
important, and Ho proposed a different approach based on the exact solution of Eq. ([1]). 
Define dimensionless bosonic operators = (x + ipx) / and similarly for Uy 

QQ- The 

one-body Hamiltonian in Eq. ([1]) is readily diagonalized with the following operators 

«± = («x T iay) and «± = ^ (4 ± ^4) (^) 

that destroy and create one quantum with frequency ujj- = l^Q and angular momentum 
±1. The Hamiltonian becomes 

Hq = 1 + u;+a^a+ + tU-oLa-, (7) 

and the angular momentum has a corresponding simple intuitive form Lz = a\_a^ — a}_a^. 

In this hehcity basis, the operator identity fliV'oo = gives the ground state (/^oo oc e"*"^/^. 
For rapid rotation — > 1, the positive- helicity normal- mode frequency uj^ = 1 — fi is small, 
whereas the negative-helicity normal-mode frequency cu- = 1 + fi ~ 2 is of order unity. The 
lowest Landau level is the set of low-lying states with with (a^a+) = m and (a^a_) = 0. 
The corresponding (normalized) low-lying excited states are 

^^{x,y) = ^=e-^-'l\ (8) 
V Trm! 

where z = x ^ iy ys, & complex variable. 



For this symmetric trap, the projector onto the lowest Landau level (LLL) follows directly 
from these normalized wave functions 

oo 

no(a;, y; x', y') = ^ Vm(a;, y) ^^{x', y'), (9) 

m=0 

and a straightforward calculation yields 

no{x, y;x',y') = ^ exp {zz'*) exp (^-^kT - ^I^T^ ■ i^^) 
The projection of any reasonable variational trial state tpv{x,y) onto the LLL basis 

Ao{x,y)= / d^r'Uo{x,y;x'y')ij^{x',y') (11) 



has the form iljyQ{x,y) = f{z) exp{—\z\'^/2), where f{z) is an analytic function of z = x + iy, 
and the remaining factor is effectively the ground-state wave function. Unless ipy already is 
in the LLL manifold, ^^^,0 typically has a reduced normalization. 

Assume a normalized linear combination of LLL solutions with ipLLL = J2m'^"i'ikm- The 
general properties of these LLL eigenstates yield the important and simple result [2] 

{L.) = {r') - 1, (12) 

where (■ ■ ■ ) means an average with iJj tj, t, - For these LLL states, the mean-field GP energy 



functional ([2]) has a simple form [3|, 

E'[,ljLLL]=n + J d\ 



22| 



(1-1]) r^\^LLL? + \g2dN\lPLLL\ 



(13) 



For numerical estimates, I use as/ Z 'o^ 10^"^ and = 5 x 10^, so that g2dN = AirNas/Z ^ 
300, in qualitative agreement with one relevant experimental coupling constant AnagN ~ 
130 {23!, where N is the number of particles per unit axial length. 

The approximate LLL ground state follows by minimizing Eq. flT3|) for fixed particle 
number and gives a parabolic density profile in the radial direction 



2 r 
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riLLLir) = \^LLL{r)\' = — ^ ( 1 - ) (14) 



„n2 \^ tj2 



Here, the mean-square condensate radius 



diverges as 1 — f2 becomes small [21]. I follow Aftalion, Blanc, and Lerner [10| and write 
1 — = for this symmetric trap. Alternatively, 1 — ^ and Eq. f[T5l) can be 

rewritten 

„o 4 In as 

which is precisely the same as Eq. (j4j) for the mean-field TP regime (as seen below, the 
situation is slightly different for a non-symmetric rapidly rotating trap). 

The validity of the LLL regime requires that the mean interaction energy Ng2dn{0) should 
be small compared to the energy gap cu- ~ 2 from the next Landau level 

= 1^ < u;- ^ 2. (17) 

A combination of Eqs. (fT6|) and f|T7|) gives the criterion to achieve the LLL regime in a 



symmetric trap 



m 



for = 5 X 10^. Equivalently, the soft frequency = 1 — Q has the restriction 1 — 17 ^ 
~ 2 X 10^^. Recent JILA experiments with larger A^ ^ 5 x 10^ have achieved 
1 — ^ 5 X 10"'^, exceeding the qualitative rotation speed for the validity of the LLL 
regime 2^. 

What is the ultimate fate of such a rapidly rotating symmetric condensate? Several 



authors 
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Q, Q, Q, Q have proposed that the superfluid LLL mean-field state 
undergoes a quantum phase transition (QPT) to a non-superfluid correlated many-body 
state. This transition occurs at a critical value of the "filling fraction" N/N^ = fc, where A^ 
is the total number of particles, A^,; is the number of vortices, and fc is variously determined 
to be between 6 and 10 (for definiteness, I take fc ~ 10). These numerical studies have 
used toroidal or spherical geometries with essentially no boundaries. To apply these ideas 
in the present case of a bounded circular condensate, note that each vortex occupies a 
dimensionless area vr for 1, so that ^ Hence the criterion for the QPT 

becomes N/R\j^j^ ^ /c, which yields the explicit critical value 



m 



4^16/2^^1.6x10-1 (19) 
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III. NON-SYMMETRIC ROTATING HARMONIC TRAP 



Now consider the more general case of a rotating non-symmetric harmonic trap 
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0, y, y, 



25| with distinct trap frequencies < Uy. In the hmit of rapid rotation, the condensate 
again becomes effectively two dimensional, and it is convenient to use the variable uj'j_ = 
(c<j^ + ojy)/2 to define both the characteristic frequency uj_[_ and the characteristic oscillator 



length ^Jh/ {Muj±). With these dimensionless units, I rewrite the squared trap frequencies 
as ul = 1 — u"^ and cu^ = 1 + z/^, where u"^ = {ujy — uj1)/2 characterizes the splitting of the 
squared trap frequencies (I follow the notation of [10]). 

A Bose-Einstein condensate in such a rapidly rotating non-symmetric trap tends to ex- 
pand along the weak trap direction oj^, and the condensate elongation apparently diverges 
as the rotation frequency Vt tends to uj^- Thus it is convenient to generalize the small pa- 
rameter as follows: = uol. — = 1 — v"^ — D?. For small e, this relation simplifies to 



A. Mean-field Thomas-Fermi regime 

It is straightforward to generalize the discussion of Sec. II. A to include a non-symmetric 
harmonic trap. The mean- field TF density remains parabolic in each direction 

\^ljTF{x,y)? = n{0) ( 1 - ^ - ^ ) , (20) 

\ ^TFx ^TFy J 

where n(0) = ji/{Ng2d) = 2/ {jiRTFxRTFy)- A detailed analysis yields the explicit expres- 
sions for the two squared TF condensate radii 

p2 , /iV^ (1 + z.^ - n^f" . /iV^(e^+2z^ 

= 'S^^Z-TZ^^ = e3/^ ' (21) 

- ^ V Z (1 + ^2 _ ^2)3/4 - 4 V Z (^2 + ■ ^ ^ 

The larger TF dimension R^fx grows continuously with increasing Vt^ 1 — z/^. In contrast, 
the smaller TF dimension R^^py grows only for Vt^ < 1 — 2z/^ (equivalently > z/^). For Vt^ > 
1 — 2z/^ (equivalently, < z/^), the smaller TF condensate dimension decreases continuously. 
As discussed below, the rapidly rotating non-symmetric system eventually reaches the mean- 
field lowest Landau level regime, requiring a different approach. 



B. Classical rotating non-symmetric trap 



With these dimensionless units, the one-body Hamiltonian now has the form 

Ho = l {pI +pI) + 1 + - ^ (xpy - yp,) (23) 

where the last term couples the x and y motion, The classical dynamical equations readily 
yield two normal modes with elliptical polarization. The two frequencies are 

u;l = l + n^Ta, (24) 

where a = + 4f2^. The plus mode with "small" frequency uj+ has positive helicity with 
polarization x + if^^y, and the minus mode with "large" frequency has negative helicity 
with polarization 2/?„x + y, where give detailed expressions for (3±. 

In the present context, the most interesting question is the behavior for fixed finite u as 
Q ujx, when = c<j^ — Q"^ is small. In this limit, both uj+ and /5+ become small, with 



1 ^ e /l — z/2 
u+ ^ and P+^-\ -. (25) 



In contrast, u;_ and /?_ remain finite, with 



^ V2V2 - 1/2 and (3^^V2\ll- — ^. (26) 

2 — 1/2 



C. Quantum states 

The bosonic quantum operators A± and A^. for a non-symmetric rotating trap correspond 
to those in Eq. ([6]) for a symmetric trap. They obey the familiar bosonic commutation 
relations = 1 with all other commutators vanishing. Nevertheless, they are rather 

complicated, and the details are relegated to Appendix A [see Eqs. (IA15P and (1A16P ]. The 
normalized ground-state wave function ipoQ{x,y) can be defined by the relations A±{pQo = 0; 



it has the explicit form 



where 



/ \ 1 ( x^ ?/2 . ac — z/2\ 



2 _ ^^(l + /3+/3-) ,2 _ + ...^ 

— — ' — — ' ^ ^ 



and 

c - (29) 

Note that c — > for a symmetric trap because in that case /5+ = /5_ = 1. In contrast, c —>■ 1 
for a non-symmetric trap as e ^ 0, since /5+ oc e — 0. 

For rapid rotation (e <^ 1), the frequencies u± have very different magnitude, with a;+ 
proportional to e, and uj_ approaching a constant of order 1. Thus the analog of the lowest 
Landau level in the present case follows with the prescription 

i^m{x,y) = iPmo{x,y) = ipoo{x,y), (30) 



where A\, has the explicit form given in Eq. flAlSI) . A detailed analysis 9| eventually yields 
the LLL wave functions 



^«.(x,y) = ^y H-[i^TTp;p:)^oo{x,y), (31) 

where Hm is the usual Hermite polynomial and 

z = x + i(3-y (32) 

is a "stretched" complex variable. It is notable that, apart from the ground state (foo, this 
wave function is an analytic function of the stretched complex variable z jsl, . Its structure 
is very similar to Eq. for a symmetric trap, and Eq. fl3ip reduces to this elementary form 
when u 0. 

The projector onto the LLL states follows directly as in Eq. (Q. Remarkably, this ex- 
pression can be evaluated analytically using a standard formula for the sum of products of 



Hermite polynomials 26|] 



Z2 2-^1 ^ =W{u,v,c) = -^ 



m=0 



C 



exp 



2uvc — {u^ + f ^) 



(33) 



where 



a[3+ z a(3+ z'* 

and z'* = x' — iB-y' . Here I focus on the behavior for fixed z/ ^ as e becomes small. 

n 

Reference [27[ uses a similar approach to study two limiting cases: a circular trap with 
z/ = 0, and a non-symmetric trap (fixed z/ 7^ 0) in the limit = uo^ (namely e = 0), when 
the condensate becomes infinitely elongated. 



When combined with the exphcit ground-state wave functions from Eq. (127|) . a detailed 
calculation yields the expression 

01 ( cazz'* \ 

no(x, y; x'y') = x{x, y) X*{x\ y') exp ( 1 , (35) 

where 

ax^ aP^y"^ . u'^xy 



Apart from the phase, this result is equivalent to that given in following a result of 
Bargmann f^]. In particular, for a non-symmetric rotating harmonic trap, the projection of 
any trial state il)v{x, y) has the form y) = f{z) x{x, y) [lot as seen from the analogous 

Eq. (HI]). 



D. Variational trial state and physical consequences 



Assume a normalized linear combination of LLL states with ipuL = J2m i^m and 



form 



9| 



1. The corresponding total energy functional in the rotating frame has the 



+ d'r 



P+X + —]WlLl\ H ^ I^LLLl 



(37) 



2^ V (3+1 ' Z 

Minimization of Eq. (!37|) with respect to IV'lllP at fixed total number yields an 



anisotropic parabolic density 

\iJLLL{x,y)\^ 

with 



, 4 NasQ , 9 ^ NasQ 



(38) 



(39) 



In the limit of rapid rotation when e ^ and a — > 2 — z/^, the squared mean-field radius 
along the weak trap direction 



Rl^ 



4^ f 2~p^ V^ IN a. 



(40) 
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diverges strongly, proportional to e In contrast, the squared mean-field radius along the 
strong trap direction 



shrinks slowly, proportional to ^/e. Both of these e dependences are the same as the leading 
factors in Eqs. fl2T]) and fl22l) describing the mean-field TF regime, but the overall numerical 
coefficients are somewhat different. 

For a non-symmetric trap with finite u, the validity of the LLL regime requires that 

Ng^MO) = 4^ ~ ^ V2V2^. (42) 

ZKxJMj 

A combination of the previous results yields 

as the criterion for the LLL regime. This expression scales as N'"^, in contrast to the A^~^ 
dependence for a symmetric trap in Eq. ( ITSl) . With a typical anisotropy parameter u = 0.5, 
the squared trap frequencies are = 0.75 and = 1.25, leading to e|j^j^ ^ 2 x 10^'^ for 
= 5 X 10^ and Z/a^ = 200. These numbers make the LLL regime seem less accessible for 
a non-symmetric trap, but reducing can increase these estimates significantly. 



IV. QUANTUM PHASE TRANSITION TO CORRELATED STATE 

As the non-symmetric trap rotates more rapidly (with = u"^ — Q^ ^ 1), the mean- field 
condensate expands in the weakly confined direction and contracts in the tightly confined 
direction [see Eqs. (HO]) and (HTI) ]. 



A. Angular-momentum criterion 

Eventually, the angular momentum arises in part from the increased moment of in- 
ertia that reflects the elongation j^, along with the more familiar vortices that dominate 
the angular momentum for a symmetric trap. This phenomenon occurs even for classical 
incompressible fluids in a rotating elliptical container [3, |2^ . Indeed, in one scenario, the 



rapidly rotating non-symmetric condensate has essentially no visible vortices lOj. 
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This situation requires a generalization of the usual filling fraction / (usually defined as 
the ratio N/N^, which is equivalent to N/R^ for a symmetric trap in the mean-field LLL 
limit). It is easy to verify that the dimensionless angular momentum per particle in a rapidly 
rotating symmetric trap is Rl/3. Thus, I here assume that the corresponding relevant ratio 
in a non-symmetric trap is / = N/{3Lz). This ratio N/L^ plays a special role in analytical 



and numerical studies of rotating bosons with relatively small |3d . l3l[ |. 

In the mean-field LLL picture with many vortices, the angular momentum per particle 
has a simple form 9| 

L. -^y. V^) + ^ (/J- - ^) (.^ - V^) - i (^- + ^) . (44) 

where (■■■) denotes an expectation value with the trial state ipiLi- It is not difficult to 
evaluate this quantity with the mean- field wave function in ( l38i) . and a detailed calculation 
yields 



12 ^ 12 y 2 V P- 
for e — > 0, where a ~ 2 — z/^. The first term grows proportional to e~^/^, the last term 
remains constant, and the second term shrinks proportional to y/e. Hence the first term 
dominates as e — and yields the approximate value 



(46) 

Compare this e~^/^ dependence with that for a symmetric trap where -Rq/^ = 



|e~^^/NaJZ from Eq. (HI). 

With this generalized definition of the filling fraction / = N/ (SL^), it is natural to suggest 
that the quantum phase transition from a coherent superfiuid state with a macroscopic 



condensate to a highly corre 



as for a symmetric trap 11, 



ated many-body state occurs at the same critical value fc ~ 10 



12 
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14| . This criterion leads to the critical rotation rate 



4 - McY^' (1 - '^Y' (2 - -'f' (^) , (47) 

which now scales as A^~^/^ instead of A^~^ for a symmetric trap [see Eq. f|T9l) ]. Take u = 0.5, 
Z/as ^ 200, ^ 5 X 10^ (as before), and fc ~ 10, which gives 



^ 3.7 X 10~l (48) 
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Note that this "angular-momentum" value of eg is somewhat larger than that for the va- 
lidity of the LLL approximation in Eq. fj43l) . Different choices for the relevant parameters 
(especially smaller N) can change these values considerably. 



B. Fragmentation criterion 

An alternative and quite general criterion for the quantum phase transition to a correlated 
many-body state relies on the following comparison 17| . Imagine the N particles distributed 
separately over the lowest single particle states in the lowest Landau level (they have 
energies juj^, where j = 0, ■ ■ ■ , — 1). This fragmented correlated state has an energy per 
particle of order Nuj^/2 and is not superfiuid. It also has small interaction energy because 
of the fragmentation. Compare this energy with the interaction energy in the lowest Landau 
level superfiuid state Ng2dn{0). The correlated many-body state should be favored when 
the fragmented energy becomes smaller than the interaction energy: ^ Ng2dn{0). 

For a symmetric trap, the small frequency is 0;+ = 1 — 1] ~ In this case, ^(0) = 

2/(7r_R|^^), and a combination with Eq. ([TB]) yields the approximate critical rotation speed 
for the quantum phase transition in a symmetric trap 

4^64^. (49) 

This value is essentially the same as that in Eq. f|T9|) based on numerical studies of the filling 
fraction. 

The situation in a non-symmetric trap is somewhat different, because Eq. ( 1251) shows 
that the small frequency in now of order e (instead of e^). Furthermore, the corresponding 
central density is n(0) = 2/{TTRxRy), so that the interaction energy becomes 

,2\ 1/4 

The correlated state has a lower energy if e < eg, where 

256 (2 -z^ / a, x2 

^ [nzJ ■ ^^^^ 



Ng^MO) ^ 2V^e { ^ ) J (50) 



4- 

Note that this N~'^ dependence is similar to Eq. (H3l) for the validity of the LLL regime 
in a non-symmetric trap. This functional dependence should be compared with that for a 
symmetric trap given in Eq. (fT9|) . With the same numbers as before, this "fragmentation" 
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criterion yields the unattainable value eg ~ 7 x 10~^. Dalibard 17| points out that this 
procedure comparing the energy of a fragmented correlated state with the interaction energy 
of a corresponding condensed state applies for other quantum critical points beyond the 
present case of a rapidly rotating dilute Bose-Einstein gas [4]. 

Which of these two criteria is correct? For the symmetric trap, the mean interaction 
energy involves the squared radius R\ii- For a non-symmetric trap, this quantity can be 
generalized either as |(-R^ + -R^), which applies for the "angular-momentum" criterion, or 
as RxRyi which applies for the "fragmentation criterion." Each gives the same approxi- 
mate criterion for a condensate in a symmetric trap. As seen above, however, they yield 
qualitatively different values for the appearance of the correlated non-superfluid state in a 
non-symmetric trap. Ultimately, this becomes an experimental question. 



V. FORMATION OF AN EFFECTIVELY ONE-DIMENSIONAL CONDENSATE 

For a non-symmetric rotating harmonic trap with uJ^ < uJ^, the mean-field LLL picture 
from Eqs. (1551) and fl5^ predicts that i?^ vanishes slowly (oc y/t) as ujI — Vt^ = ^ 0. 
Ultimately, such behavior must cross over to a different dependence, and Aftalion et al. flo| 
suggest the formation of a quasi one-dimensional regime with no visible vortices. This idea 
arises from the LLL projector in Eqs. fl35|) and fl36l) . where the factor xi^^v) implies a 
Gaussian dependence on the tightly confined variable y. 

If the trial function is a narrow Gaussian ipv{y) oc exp(— ?/^/2cr^) with mean width a, the 
projected function tp^o has minimum energy E' when a — > 0, yielding a projected function 

i^voiy) oc exp{-y^/2yl^-J, (52) 

with y^jjj ^ Q/a(3- [l^. Hence it is natural to assume a crossover from the mean-field 
parabolic state to a quasi one- dimensional state when ^ Cy^^^, where the choice of the 
numerical factor C is somewhat arbitrary. 

For rapid rotation with e — 0, the characteristic squared dimension ?/^-^ approaches a 
constant value 

= ^ VWT^' ^^^^ 

Assuming that the mean-field parabolic density profile in Eq. fl38p can indeed occur in the 
lowest Landau level, the crossover from the parabolic parameter Rl to the Gaussian i/^^ 
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occurs at a corresponding rotation rate 



,2 ^^i2-z^,^. 



1024 d - z/2)3 V A^a 



With C = 10, my previous numerical parameters {u = 0.5, = 5 x 10^, and Z/a^ = 200) 
imply e^j.ogs ^ 1.0 x 10~^, which might well be attainable. Recent numerical studies of the 
formation of a single vortex row suggests that the condensate's order parameter remains 
Gaussian in this situation, so that the constant C may well be as large as 20 32|]. In this 
case the crossover parameter becomes Ccross ~ 1-6 x 10^^, which should be readily observable. 

Figure 2 of Aftalion et al. 1^ displays an effectively one-dimensional condensate with no 
visible vortices, but they assume a much smaller value Nag/az ~ 3, where is the extension 
of the wave function in the z direction 33|]. Such a weak coupling constant would require 
a smaller particle number N, since is ~ a few nm and a typical quasi one-dimensional 
condensate has a transverse dimension ^ 1 /xm. 



VI. CONCLUSIONS AND OUTLOOK 



A rapidly rotating symmetric two-dimensional condensate expands radially and has a 
large squared radius -RI^^l proportional to 1/e, where = 1 — Visible quantized vortices 
provide the total angular momentum L^. In oscillator units, each vortex occupies an area 
~ TT, so that the number of vortices is A^^ ~ "^^^ criterion for a quantum phase 

transition to a correlated state is usually given as N/Ny f^. Taking fc ~ 10, I find 
that the critical rotation rate for this quantum phase transition is ~ 2 x 10^'^, which is 
experimentally challenging but perhaps feasible. 

The LLL wave function for a symmetric trap ensures that (L^) = (r^) — 1, and the 
parabolic density profile from Eq. (fT^ yields the approximate result (L^) ^ \R\ll 
large R\ll- Hence the criterion for the quantum phase transition for a symmetric trap is 
equivalent to N/Lz ~ 3/c, which is more general than N/N^ in cases where the angular 
momentum arises from other mechanisms (such as in a non-symmetric trap — see below). 

A non-symmetric trap has ujy — uj'^ = > 0, and it is now natural to define the small 
parameter as = — Q"^. For rapid rotation, the condensate expands along the weak x 
direction and contracts along the tight y direction. In one scenario for a quantum phase 
transition, the resulting asymmetry induces a large moment of inertia and a correspondingly 
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large angular momentum proportional to e~^/^. Note that this dependence is more 
singular than that for a symmetric trap, where oc e~^. As a result, the relevant ratio 
N/ Lz for a non-symmetric trap decreases more rapidly than for a symmetric trap, favoring 
the quantum phase transition to a correlated state. 

A second scenario for the quantum phase transition involves a fragmented state |171]. In 
contrast to the angular-mementum criterion, this picture implies that the quantum phase 
transition is more difficult to achieve for a non-symmetric trap. 

Instead of a quantum phase transition, a third scenario predicts a smooth crossover to a 
quasi one-dimensional condensate with no visible vortices 



of symmetric traps 



11 



12, 



13 



lOj . Previous numerical studies 



14| relied on geometries like tori or spheres that have no 



boundaries. More general numerical studies that include the trap anisotropy would be 
highly desirable. 



351,1361,1371. 



One intriguing experimental possibility is to use synthetic gauge potentials [34 , 
Such light-induced effective vector potentials can mimic a rotating system while remaining 
in the laboratory frame. It would be very interesting to have such experiments on a non- 
symmetric trap, which could help decide among these various scenarios. 
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APPENDIX A: CANONICAL TRANSFORMATION FOR ROTATING NON- 
SYMMETRIC TRAP 

The Hamiltonian in Eq. fl23l) is quadratic in the coordinates and momenta, so that a 
linear canonical transformation can diagonalize it, but the details become quite intricate. 
Introduce a four-component vector v'^ with elements x, ?/, p^i Py, where T denotes a transpose 
of a vector. The one-body Hamiltonian in (|23l) then has the form Hq = ^v^Tiv. In the 
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present case, the symmetric 4x4 matrix Ti is 



H 





^] 1 

-1] 1 



V 



(Al) 



The Hamiltonian equations of motion are equivalent to 

dH 



V = (T^— = (tHv, 

ov 



(A2) 



where the 4x4 matrix 



a 



(A3) 



/ 1 o\ 
1 
-10 
\ -10 0/ 

enforces the symplectic structure of Hamilton's equations. I seek normal modes with time 
dependence oc e"""^*. In this case, the dynamical equations (]A2p become a linear matrix 
equation 

Tiv = uiav, (A4) 

where ia is a Hermitian matrix with (ia)^ = ia. 

The eigenvectors Uj obey the corresponding eigenvalue equation 



TCuj = ojjiauj 



(A5) 



with two independent real roots uo± given in Eq. ( 12^ . The resulting eigenvectors are or- 
thogonal with respect to the matrix ia that serves as a metric 



u^jiauk = 5jk. 

Straightforward but lengthy calculations yield the normalized eigenvectors 

/ 1 \ / ,7, \ 



(A6) 



and u_ 



A_ 



1 



(A7) 
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where A_|_^ = 1 + l3+l3-UJzf/uj±. 

Use these eigenvectors to construct the 4x4 matrix 

/ Ml U U*. U* \ 

• (A8) 

V i i i i / 

The orthogonahty relations (]A6I) are equivalent to the single matrix equation 

M^iaM = T, (A9) 

where r is a 4 x 4-dimensional generalization of the familiar Pauli matrix cr^ with diagonal 
elements (1, 1, —1, —1). In addition, the matrix J\f satisfies the relation [compare Eq. flA5l) ] 



KM = iaAfrO, (AlO) 

where O is a diagonal matrix with positive elements uj+, uj-,uj+,uj-. The normal coordinates 
A_ can be defined through the linear transformation 

V = AAA, (All) 

where the 4-dimensional vector A has the explicit form (here given as the transpose) 

A^= (A12) 

This transformation also diagonalizes the original Hamiltonian 

H = ^A^OA = ^uj+ (^AIa+ + A+Al^ + ^cj_ (^aIa_ + A_aI^ . (A13) 

Equations (1A9P and (lAlip yield the explicit form of the normal coordinates 

A = rAfhav, (A14) 

in terms of the eigenvectors u± and the original canonical coordinates x,y,px,Py A detailed 
analysis yields the rather complicated expressions 

A+ = [(^+ + np+)x + p+py - i{n + p+u+)y + ipx] , (A15) 

A- = ^= [(cj_ - np^)y + (3-Px + i{n - p_uj.)x + ipy] , (A16) 
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that generalize the results ([6]) for a rotating symmetric trap to the case of a rotating non- 
symmetric trap. In the symmetric limit (z/ — > 0), a+, but —ia^ = {ay — 
iax)/\/2, which is merely an additional canonical transformation. 
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